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We present predictions for the Z-boson pt-spectrum at Tevatron within the framework 
of unintegrated distributions evolved according to evolution equations recently proposed 
by us. We discuss the dependence of the results on the choice of non-perturbative pa- 
rameters, the coupling constant and the impact of soft gluon resummation. 
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1. Introduction 

Transverse momentum dependent, or equivalently, unintegrated distributions are 
currently object of intense research activity. The motivation for such an interest re- 
Ues on their wide range of apphcabihty, from spin physics to jet observables in high 
energy colUsions at LHC. Their correct formahzation in quantum field theory, in 
the present case in quantum chromodynamics, is however far from being triviaW In 
order to investigate various properties of unintegrated distributions, detailed calcu- 
lations in various gauges have been performed— - '. The structure of additional rapid- 
ity divergences is understood and ad hoc subtraction scheme has been proposeci^. 
Their factorization properties in hard processes have been also investigated and a 
factorization theorem has been given^. A point which, to date, escapes a rigorous 
answer is how these distributions behave as long as the scale which characterizes the 
hard process is varied. If this behaviour were known we would be able to relate to 
each other results coming from different experiments, possibly at different energies. 
The latter possibility is therefore of great phenomenological importance. Although 
a definitive answer to this question, especially in the light of new developments in 
the field, is absent in the literature, there have been however some attempts. In 
particular evolution equations for unintegrated distributions were proposed in the 
unpolarized time-like casJ^ and very recently extended to space-like kinematicJ^. 
In a subsequent phemonelogical studj^, performed in the context of semi-inclusive 
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deep inelastic scattering, it was shown that a reasonable description of data could 
be obtained once unintegrated evolution equations were solved with suitable, but 
motivated, initial conditions and assuming factorization for the cross-sections of in- 
terest. This result has stimulated us to apply the same formalism to Drell-Yan type 
process in hadronic colhsions. The pt-spectrum of the gauge boson has, in fact, a 
rich structure and manifests many perturbative and non-perturbative features of 
the underlying theory. In particular techniques for the resummation of the pertur- 
bative series in the multiple soft gluon emission limit were first developed for this 
prototype 

observabltPllnl. 
2. Unpolarized evolution 

We briefly summarize the basic ingredients of fct-evolution equations. Let us consider 
parton emissions off a active, space-like, parton line in ladder approximation. In the 
coUincar limit, at each branching, the active parton increases its virtuality and 
acquires a small relative transverse momentum with respect to the parent. These 
iterated emissions generate therefore an appreciable transverse momentum, up to 
the order of the hard scale in the process, which adds to the non-perturbative 
one due to Fermi motion of the parton in the parent hadron. CoUinear emissions 
give however leading logarithmic corrections to cross-sections when the transverse 
momenta are ordered along the ladder and can be resummed to all orders by using 
DGLAP evolution equationJll] ^j-^g unintegrated casc"^*^, unlike DGLAP case, 
the integration on relative transverse momenta at each branching are left undone 



as indicated by the additional (Pl± integration. In the splitting k k + p depicted 
in Fig. Il]), the parton k carries a fraction u of the fractional momentum of k. With 
this notations the following mass-invariant constraint can be derivered: 



If one assumes that the virtualities increase along the ladder, fc^ 3> fc^, and on-shell 
partons are emitted, = 0, the last two terms in eq. ([2]) can be disregarded. In 
these limits, setting —k'^ = one obtains = (1 — u)Q^, which can be found in 
the i5-function in eq. ([T]). The transverse arguments of Tp on r.h.s. of eq. ^ are 
derivered by taking into account the Lorentz boost of transverse r nom enta from the 
emitting parton k reference frame to the interacting k parton on4^. In particular, 
the transverse momentum k± of the parton which undergoes the splitting can be 
expressed as follows 



Q 




(1 - u)k^ + u{l - u)k'^ - up^ . 



(2) 



fe^ = {k±_ - li_)/u. 



(3) 



Trento proceedings 3 




Fig. 1. The incoming proton on the bottom left is assigned to have unitary momentum and 
defines the reference axis respect to which the transverse momentum fc^ and fc^ are defined. 
The small blob represents the iteration of emissions in the parton ladder of whose only the last, 
fc — > fc + p, is explicitely shown. lj_ is the relative transverse momentum between parton k and 
p. The interacting parton k enters the hard scattering vertex indicated by the blob on top of the 
diagram. 

Unintegrated parton distribution functions J-p {xb,Q^ , k±_ ) in eq. ([1]) give the prob- 
ability to find, at a given scale , a parton i with longitudinal momentum fraction 
Xb and transverse momentum fej_ relative to the parent hadron, see Fig. ([T]). Pji{u) 
are the space-like splitting functions. The unintegrated distributions fulfil the nor- 
malization condition: 

J £k^rp{xB,Q\ fe±) - fp{^B,Q^) , (4) 

where fp are ordinary parton distributions. It is important to remark that, given 
eq. (Ill), performing the (Pk± integration on both side of eq. ([T]), we recover the 
integrated evolution equation^iiJ for fp. 

3. The pt-spectrum of Z boson at Tevatron 

The analysis described in the following is inspired to a similar on^^^, in which how- 
ever a different kind of evolution for of unintegrated distributions is assumed. The 
differential cross-sections for Z-production at rapidity y and transverse momentum 
Pt is given by 



.F,(a;l,^^fci4)J^g(x2,M^fc±,2) + (l<->2)J . (5) 

This formula can be understood as the counterpart of the standard factorization 
formula for Drell-Yan type processeJI^. The weak charges are denoted by Wq, the 
parton momentum fractions are evaluated through xi^2 = mT / \/s exp(±y) with 
rriT = \/ -I- being the transverse mass of the gauge boson. The hadronic 
collision energy is ^/i ~ 1.8 TeV and ctq = ^'^3/'™^ • Parton distributions are 
evaluated at the scale /i^ — m\. The differential cross-section, eq. ([5]), is then 
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Z boson p spectrum at Tevatron 
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Fig. 2. Predictions for the Z-boson pt-spectrum along with Tevatron data. 



integrated according to the experimental cuts-^ and accounted for the Z branching 
ratio into electrons, BR{e^e~) = 0.033632. The initial conditions for fct-evolution 
equation at the min imum scale Qq — 5 GeV^^ are chosen as a product of longitudinal 
parton densitieJl^ times a gaussian transverse factor with x- independent and flavour 
independent widtlPl 



J^p{xB,Ql,k±) ^ fp{xB,Ql) 2 — r ^ = 1^1, 9- (6) 

We evolve light flavours onljC^I and effects of heavy flavours are included in the 
running of the strong coupling evaluated in leading logarithmic approximation. We 
tune the parameters appearing in the intial conditions, eq. ^ , as well as the strong 
coupling at the Z-boson mass, as{m\), to data. During this procedure we have 
observed two peculiar features: setting the quark intrinsic momentum {k\^ ^) to 
larger values shifts the position of the maximum towards higher pt while lowering the 
value of the coupling, as{M'^), do overestimate its height and viceversa. The tuning 
procedure gives for the quarks (and antiquarks) an intrinsic momentum (fc^ q)~'^ 
GeV^. Quite interestingly the same amount of intrinsic transverse momentum is 
also required by Monte Carlo programi^^l when used to predict the same process. 
We will discuss this large value, well above the one expected by Fermi motion, 
after the discussion of soft gluon resummed results. Predictions are, as expected, 
almost insensitive to the gluon intrinsic transverse momentum so that we fix it at 
i^i. g)~^ GeV^. In order to have a reasonable agreement with data, a large value 
of the coupling is required, namely as{M'^) = 0.150. The high value for Q;s(Mf ) 
indicates that a more rapid evolution is necessary, especially in the low pt region. In 
order to have a better description of the pj-distribution tail, we use the full invariant 
mass-constraint at the branching vertex, eq. ([2]) , still however considering emission 
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Fig. 3. Soft gluon resummed predictions for the Z-boson pt-spectrum along with Tevatron data. 

of massless partons, = 0. The predictions within these settings are shown, along 
with experimental data^^, in Fig. ([2]). The main features already observed in the 
analysis^ of semi-inclusive deep inelastic scattering appear here again: the use of 
/ct-dependent distributions allows one to describe gauge boson p^-spectrum without 
resorting to any artificial procedurJI^l to match higher order pQCD corrections at 
large pt and non-perturbative predictions at small pt, possibly corrected for soft 
gluon emissions. This quite interesting feature is the result of taking in full account 
the transverse kinematics in the proposed evolution equations. As is well known, 
the low Pt part of the spectrum is sensitive not only to pure n on-pert urbative effects 
but also to effects coming from multiple soft gluon emissions.^^Min! jn this case the 
fct-evolution equations can be slightly modified in order to resum logarithms of soft 
naturJI^. In the non-singlet channel the resummed evolution equation reads 



Q 



^drp{xB,Q^,k±_) _ f durasiQ^il-u)) 



5Q2 



J ^s[il~u)Q^-ll)Tj.[^,Q 



2 



In the previous equation Pqq{u) denotes the unregularized splitting function. The 
resummation of leading soft logarithms is performed by changing the argument of 
the running coupling from the virtuality to the relative transverse momentum of 
partons at each branching, as(Q^) aJl'^i )■ We wish to note that eq. ^ can be 
easily recast in impact parameter space^. Taking advantage of the solution of the 
/ct-evolution equations in the soft limilP^, is then possible to recover the well known 
results for the perturabtive form factor in double logarithmic approximatiorP. The 
off-diagonal term in the Altarelli-Parisi splitting matrix are assigned to have the 
standard coupling, as{Q^), since no soft enhancement is present in the q —> g{g) 
and g — s- q{q). The soft gluon resummation can be also performed in the gluon 
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Z boson p spectrum at Tevatron 
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Fig. 4. Effects on the predicted distribution by the variation of the quark intrinsic momentum. 
Curves are shown for the following values : {k^ ^)=2,3,4 GeV^. 

channe i^ an d the whole formalism extended up to next-to-leading logarithmic a,c- 
curacy D^^l'^^l . The latter improvements however are not still implemented. Despite 
for u — !■ 1 virtual and real contributions exactly cancel, as guaranteed by the plus 
prescription, the rescaled coupling can be sensitive to how the infrared limit is ap- 
proached. We adopt here the simplest modeP, e.g. a freezed couphng as{k^ + {g^}) 
being g^ = 0.5 GeT^^ the freezing scale. We stress that more refined prescription 
for the infrared behaviour of the strong coupling could be usecfl^. In Fig. Q we 
show the predictions which include soft gluon resummation as we have described 
above. The tuning procedure gives an intrinsic momentum (fc^ g)=3 GeV^, and 
a much lower value for the coupling, as(M^) = 0.120. Moreover the insensitivity 
to gluon parameters persists so that we still fix the gluon intrinsic momemtum to 
(fc^ g)=l GeV^. The inclusion of higher order terms in the perturbative calcula- 
tion provided by the resummation has therefore the effect of strongly reduce the 
coupling and slightly reduce the quark intrinsic transverse momentum, not still in 
the range of what expected by Fermi motion. In the resummed case, the sensitiv- 
ity to the quark intrinsic momentum is shown in Fig. (j4]). The fact that a large 
amount of intrinsic momentum is required in the description of the Z-boson data 
is known in the literature and it persists even when soft resummation is pushed to 
next-to-leading logarithmic accuracjEOl There are however attempts to solve this 
probl em, either using x-dependent (i.e. energy dependent) non-perturbative form 
factor^ or introducing modifications to the coupling constant The former ap- 
proach is particularly suitable when combined analysis of low energy Drell-Yan and 
Z 01 W data are performed. A possible energy dependence of the non-perturbative 
form factor used in that works, and here embodied in the transverse part of the 
initial conditions, eq. ([6|), is therefore of great importance in view of LHC Drell-Yan 
physics program. In order to asses the reliability of all these results a careful anal- 



Trento proceedings 7 



ysis of the uncertainty due to scales and parameters variation is presently under 
way. Furthermore an analysis of low energy Drell-Yan is also planned as well as the 
implementation of soft gluon resummation in the gluon channel. 

Conclusions 

We have found that the predictions based on the proposed evolution equations 
for the unintegrated distributions are able to reproduce Tevatron data on the pt- 

spectrum of the Z-boson. However the imrcsummcd results arc characterized by 
large values of the quark intrinsic transverse momentum and coupling constant. 
The inclusion of soft gluon resummation strongly reduces the value of the coupling 
and brings it more close to world average. On the other hand it only slightly re- 
duces the large values of the quark intrinsic transverse momentum with respect to 
unresummed predictions, leaving therefore open the problem of the nature of the 
latter. Owing to these results and within the evolution equations scheme adopted 
in this analysis, wc conclude that non-perturbative parameters differ significantly 
from the one used in our semi-inclusive DIS data analysis. 
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